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We investigate low- lying singlet and triplet excitations and the magnetization process of quasi- ID 
spin systems composed of tetrahedral spin clusters. For a class of such models, we found various 
exact low-lying excitations; some of them are responsible for the first-order transition between two 
different ground states formed by local singlets. Moreover, we find that there are two different kinds 
of magnetization plateaus which are separated by a first-order transition. 
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I. INTRODUCTION 

Recently, there has been substantial interest in the 
properties of strongly frustrated low-dimensional quan- 
tum spin systems such as S = 1/2 Heisenberg models 
on the kagome and pyrochlore lattices. In such systems, 
various ground-state phases, e.g. valence-bond crystals, 
RVB, long-range ordered phases etc, appear as we vary 
the control parameters. In particular, the possibility of a 
spin-liquid phase with unconventional singlet excitations 
would be interesting!!!. 

In the following, we consider a simple model system 
built up from tetrahedral clusters of S = 1/2 (we label 
each tetrahedron by Latin indices j) with the Hamilto- 
nian (see Fig.|l|) 
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These tetrahcdra form a chain-like structure and interact 
with each other by the following coupling: 

"Hi = ^[^3 (Sj,2'Sj + i,3 + Sj^-Sj + i^i) 

j 

+ Ji (Sj : 2-Sj + i ; i + Sj^-Sj+i^)] . (2) 

For J2 = J3 = J4, this model reduces to a known model of 
a frustrated spin ladder with diagonal couplings ('gener- 
alized Bose-Gayen model trtil). The choice J3 = J4, ^ J 2 
(symmetric model, hereafter) introduces explicit dimer- 
ization in the leg direction. As can be seen in Fig. an 
assembly of decoupled tetrahedra described by Tto has 
singlet modes with energies much lower than the singlet- 
triplet gap. Therefore, the model with J±, J2 3> ^3,^4 
would provide a good starting point to study the uncon- 
ventional properties mentioned above. 

The study of the model described by the Hamilto- 
nian Htot = 'Ho + "H-i W a s inspired by the discovery 
of the tellurate materialsa Cu2Te 2 5 X5 (X= CI or Br) 
and the phase diagram, triplet excitations, and the opti- 
cal spectrum of this model have been investigated quite 



recentlya using both numerical diagonalization and the 
bond-operator mean-field approximation. In these mate- 
rials, Cu 2+ ions form S — 1/2 tetrahedra, which are con- 
nected with each other by Te-0 coordinations. Although 
the coupling between the tetrahedra is not so simple, the 
crystal structureQ (P4) suggests that the model Hamilto- 
nian 7itot with J3 = J4 is one of the simplest candidates 
to describe the tellurates (the chain axis is parallel to 
c-axis of the tellurates). On the basis of experimental 
results obtained far static magnetic propcrticscl and for 
Raman scatterings, it was argued that the parameters 
Ji py J 2 and J3 = J4 <C Ji may be appropriate for the 
two compounds. We will also include some results for a 
more general case with J3 =/= J4 (hence we call it gener- 
alized model in the following) which may clarify to what 
extent our results are general. 

Our aim here is to investigate (i) general excited states 
and their relation to the ground-state phase transitions 
and (ii) the magnetization plateaus. In what follows, 
we will use eigenstates of the four S — 1/2 spins on a 
single tetrahedron in the notation as shown in Fig.^. The 
ground state of the chain of non-interacting tetrahedra 
J3 = J4 = is obtained as a sequence of tetrahedra in 
state '1' for J\ > J 2 and as that of tetrahedra in state '2' 
for J\ < J 2 . 

We conclude this introduction with a few comments 
on the ground state phase diagram for the symmetric 
model (including the coupling J3 = J4) which was in- 
vestigated already in Ref. ^|. The main purpose of this 
part consists not in obtaining the phase diagram itself 
but in demonstrating how the low-lying singlets dictate 
the transition. To discuss the structure of the phase dia- 
gram from the viewpoint of the singlet spectra, we use an 
effective Hamiltonian derived by paying particular atten- 
tion to low-energy singlets. This approach gives a simple 
and clear picture of the transition between singlet phases 
and will be useful also for a discussion of the low-energy 
singlet dynamics. 

The classical model (S f 00) has two different kinds of 
antifcrromagnctic phases separated by a line J\ = J 2 + Js- 
Along the transition line, the ground state exhibits a huge 
degeneracy, which is reminiscent of what occurs to the 
classical pyrochlore antiferromagnetsEl. As was already 
pointed out in Ref. |[ the symmetric model (S = 1/2) 
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preserves the basic propertyi of the generalized Bosc- 
Gayen model: the total spins S^i + Sj 2 and S 0t s + Sj.4 
on individual dimcr bonds (Ji) are well-defined quan- 
tum numbers. This allows a simple classification of the 
eigenstates of 7itot by specifying these quantum numbers 
for all 27V (TV: the number of tetrahedra) dimer bonds. 
When all dimers are occupied by triplets, the Hamilto- 
nian 7itot reduces to an effective spin-1 chair 
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k: odd 



where denotes an effective spin-1 operator S, k+i + 
S 3 fc+i (for k: odd) or S 2 k + S 4 k (for k: even). On 
the other hand, interactions J 2 and J3 effectively vanish 
when all dimer bonds are occupied by singlets. 

It was argued in Ref. that the quantum phase dia- 
gram contains three different phases: rung-dimcr (RD), 
HaldaneEl (H) pad spin-1 dimerized (SID, or plaquctte 
singlet) phaseollJ; the first one is characterized by the 
formation of local singlets on dimer (Ji) bonds. On the 
other hand, the latter two are_phases of the model 7i 1 
(eq.(||)) and are distingjiisheaO according to whether the 
string-order parameteiO is vanishing or not. Note that 
we can use a variational argument to show rigorously 
that RD phase actually realizes at least for J\ > J 2 + J3 
(J2, J3 > 0) in the sense that no admixture of singlets '1' 
and '2' occurs. 

To investigate the effect of the coupling J3 between 
tetrahedra on low-lying singlets analytically, we derive an 
effective Hamiltonian acting on the 2 A '-dimensional sub- 
space spanned by two nearly degenerate states '1' and '2'. 
At 'site'-j (i.e. j-th tetrahedron), we define an Ising spin 
with a value +1 ( — 1) when the tetrahedron is in state 
'2' ('1'). The coupling between tetrahedra is taken into 
account by degenerate perturbation theory and the re- 
sulting effective Hamiltonian is given by a ferromagnetic 
Ising chain in an external field: 
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with 



JlM — 7T7 - 7 ^IM — J-2 — J\ + 7T7" ■ (5) 

We have different ground states according to the sign of 
the effective magnetic field hiM- For hiM < the Ising 
spins align downward and the 'rung dimer (RD)' phase 
realizes, while a positive value of hiM makes Ising spins 
point upward to form the 'spin-1 dimer (SID)' phase. 
The condition hiM = determines the line of first-order 
transition between the RD- and the SID phase. The 
mechanism of the transition will be discussed in the next 
section from the viewpoint of the excitation spectra. In 
Fig.^, we show the transition line obtained above by a 
solid line. This result implies that at the symmetric 
point J2/J1 = 1 a small perturbation J3 resolves the 



huge degeneracy and SID phase is selected as a unique, 
spin-singlet ground state. 

The explicit form of the ground state is known ana- 
lytically in the whole RD phase (all tetrahedra in state 
'1') and only on the line J 3 = 0,J 2 > Ji in the SID 
phase (all tetrahedra in state '2'). In the remaining part 
of the phase diagram, the ground state is known numer- 
ically from— the study of S = 1-chains with alternating 
exchange&L3. In particular the separation line between 
the Haldane and the SID phase is given by J3/J2 ~ 0.6, 
where a second-order transition described by the 9 = tt 
0(3) non-linear sigma modclQ occurs. The quantum crit- 
ical point, where all three phases meet and the gap of the 
first order transition disappears, is of particular interest, 
but will have to be treated beyond perturbation theory. 
In section [V, we briefly discuss the effect of interchain 
couplings in conjunction with three-dimensional order- 
ing. 

As an independent check of our method, we also de- 
termined the RD-S1D boundary by adapting numerical 
data of Ref. |l0| (open squares in Fig. ^. Note that this 
is essentially the same as that given in Ref. ||). The re- 
sulting ground-state phase diagram is shown in Fig. [| 
as a function of J2/J1 and J3/J1 with J\ as the en- 
ergy unit. The merit of this representation is that the 
phase diagram clearly exhibits the symmetry under the 
exchange of J 2 and J3. In the inset, we also show the 
same data in the (J3/J2, J\l J2)flplane, the parametriza- 
tion used in a recent literature^, in order to facilitate 
a comparison. Although we used a completely different 
method, the result is consistent with the-.one obtained by 
the bond-operator mean-field approacm. The inclusion 
of higher-order terms hardly changes the boundary and 
we may expect that the convergence of our calculation is 
good. We would like to stress here that our simple effext. 
tive Hamiltonian (Q) not only yields a fairly good resultE3 
but also gives us a clear picture of the transition as will 
be described in the next section. 

For the generalized model, the phase diagram is not 
symmetric and the relation to spin-1 chains is no longer 
useful. For small J3 (or equivalently, J4), however, the 
mapping to the Ising chain (||) goes in a similar manner 
to give 
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Setting Ji = J 2 , we obtain a negative value hiM = 
— jf/(12Ji). Contrary to the previous case (symmetric 
model), this implies that for perfect (Ji = J 2 ) tetrahe- 
dra the inter-tetrahedron coupling selects not the SID 
ground state but the RD one. The results obtained in 
this section lead us to conclude that if the value J\ = J 2 
is reliable for the real compounds CU2TC2O5X2 (X=C1, 
Br), they should be treated in the SID phase. 
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II. EXCITATIONS 

Because of the fully frustrated tetrahedral structure, 
there are various types of singlet and triplet excitations. 
Among them, singlet excitations within the triplet gap 
are of particular interestEjolIj. Main results of this sec- 
tion are summarized in Figs. |[ [| and Table I. 

A. Excitations in the rung-dimer phase 

1. Singlet excitations 

As can be easily seen from FigJ|, low-lying singlet de- 
grees of freedom do exist in the neighborhood of the point 
J\ = J2 which corresponds to perfect tetrahedra. The 
lowest singlet in the rung-dimer (RD) phase is created 
by promoting one of the RD tetrahedra to an SID sin- 
glet (state '2' in Fig.[|). For the symmetric model, this 
is an exact eigenstate with energy 2(Ji — J 2) since the 
interaction Hamiltonian Ji\ annihilates this state. For 
an estimate of the interaction effects (oc J3), the effective 
Hamiltonian (Q) can be used and leads to the results: 

App S = 2(Ji — J2) for the symmetric model, 

for the generalized model . (7) 

Alternatively this state can be viewed as a singlet bound 
state made up of two dimer triplets in the dimer region 
(J 1 ^ J2). The expression for the generalized model 
(second eq. of (Q)) is, of course, not exact, but we can 
show, at least in a perturbative sense, that this excitation 
is completely localized (i.e. dispcrsionless) also in the 
generalized model. 

What is more interesting is that there exist multipar- 
ticle bound states which is given by n(> 2) successive 
SID singlets; ferromagnetic interaction in the effective 
Hamiltonian provides the attraction between these par- 
ticles (the binding energy is — 4(n — 1)| Jim|)- For both 
symmetric- and generalized models, the total energy of 
this kind of bound states is given by 

AXund=A^ s + 2(n-l)|/nM|. (8) 

Thus, when the phase boundary is approached from the 
RD side (hiu — > —0), all these bound states collapse onto 
the lowest singlet (while a small gap of order J\ — J2 
remains between the singlet ground state and the lowest 
singlet excitation-see Fig.||). This collapse triggers phase 
separation and leads to the first-order transition from 
singlet RD to singlet SID. In particular, the largest one 
with the gap A^ bound = 2N\hiu\ hits the ground state 
at the transition point and after the transition these huge 
bound states constitute low-lying excitations of the new 
phase. Contribution of these low-lying singlets to such 



physical quantities as specific heat can be calculated by 
using the solvable Hamiltonian Wising- 

On top of them, there are several singlet bound states 
composed of two triplet tetrahedra. For example, a sin- 
glet combination of two triplet tetrahedra (state 'a' and 
'/?') connected by a weak (J3) link has an exact energy 
2Ji — 2J3, which lies between the two-triplet threshold 
2Ji and the elementary singlet A^ g . 

2. Triplet excitations 

In this subsection, we discuss several exact magnetic 
excitations for the symmetric model. Although most of 
the following results hold also for the generalized model, 
the excited states are no longer exact. 

The simplest such excitation is an immobile dimer 
triplet excitation created by replacing one of the singlet 
(state '1') tetrahedra by a tetrahedron in state 'a' or '/?'. 
The exact excitation energy is J\ , independent of J2 and 
J3. These are nothing but dimer-triplet excitations of the 
standard two-leg ladder with Bose-Gayen type couplings. 

Another triplet with energy 2Ji — J2 can be created 
by promoting one tetrahedron from state '1' to state 
'7', composed of two dimer triplets. In the dimer limit 
( Ji 3> J2, J3), this state can be viewed as a triplet bound 
state made up of two dimer triplets (with binding energy 
— J 2)- If two dimer triplets are bound on a weak (J3) 
link, then they form another exact bound state with en- 
ergy 2Ji — J3, corresponding to adjacent tetrahedra in 
states 'a' and '/?'. By a logic similar to that used in the 
previous subsection (section II.A.l), we can show that 
these triplets are completely localized (i.e. dispersion- 
less) even if we relax the condition J 3 = J4. 

B. Excitations in the S — 1 dimer phase 

1. Singlet excitations 

Starting from the ground state of the SID phase, ele- 
mentary singlet excitations for the symmetric model are 
obtained by changing some tetrahedra to state '1'. Since 
the intcrtetrahedra coupling Tii annihilates links with 
at least one tetrahedron in state '1' on their edges, this 
change results in a sequence of finite S = 1 chains. The 
excitation energy, however, cannot be calculated analyti- 
cally because of quantum fluctuation coming from S — 1 
segments. For the lowest singlet excitation, obtained for 
one tetrahedron in state '1', the energy is obtained in 
second-order perturbation as (the second-order correc- 
tion in the following expression is the contribution from 
the open ends of two S = 1 chains) 

A^g=2(J 2 - J 1 )+4J 3 2 /(3J 2 );. 

We have A g ™p > in the SID phase. It is easy to verify 
that Ag-j 1 ^ coincides with the first equation of ([?]) at the 
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transition point |/iim| = 0. 

As in the case of the RD phase, we can consider several 
multiparticle states made up of state '1' tetrahedra (both 
scattering states and bound states). Among them, the 
most important is an rt-particle bound state. The energy 
of immobile n-particle bound states can be calculated in 
perturbation theory as 

A™§ + 2(n - l)hm ■ 

Again, when h\M ~ 0, the binding energy is so large 
that all multiparticle bound states have the same energy 
Agjjj- This is the origin of instability towards the first- 
order transition from the SID to the RD. It is important 
to note that all the above excitations are not included in 
the effective spin-1 chain, but correspond to an internal 
degree of freedom on a pair of dimer spins. 

2. Triplet excitations 

Two different types of triplet excitations exist. The 
first one is essentially a single tetrahedron in state- 
'a'('/?') in the SID background; it creates a single dimer 
singlet in the sea of dimer triplets and is not included 
in the usual '5 = 1 chain'. This is highly localized be- 
cause an 'unpaired' spin-1 object appearing at the edge 
of the effective 5=1 chain (with an odd number of effec- 
tive spins) can hardly move due to strong dimerizationE3. 
The energy is given by 

A^ = -^i+2J2 + 0(jf). 

Triplet excitations of the second type are contained in 
the excited states of '5 = 1 chain' and are obtained by 
promoting a single tetrahedron to state-'7'. The excita- 
tion discussed in Rcf. | by a mean-field approximation 
etc. is of this type. Contrary to the '7' tetrahedron 
in the RD phase, this magnon excitation can propagate 
freely due to the 'background' of dimer triplets and the 
dispersion is given by 

triplet/ \ 
W type-2(9) 

/ 8J| \ /4J 3 2Jf\ 4J| 

= ( j2 + 2d) -{-t + d) cosq - d cos2q 

(9) 

Note that this has a relatively large bandwidth of 0.92 Ji 
(for J\ = J 2 and J3/J1 = 0.3). This fact shows that 
the effect of intertetrahedra couplings on low-lying exci- 
tations is drastically different according to the spin back- 
ground. This difference may be crucial also in consider- 
ing possible scenario of three-dimensional ordering (see 
section[Fv|). 

Furthermore, by solving a two-magnon problem 
cxplicitlyE3, we found singlet- and triplet bound states 
around q — tt, whose binding energies are a few percent 



of the triplet gap. The singlet one agrees qualitatively 
with that pointed out in Ref.|^ in conjunction with the 
Raman spectrum. The detail will be reported elsewhere. 

Because of the large bandwidth, a crossing between 
the type-2 mobile triplet and the lowest gapped singlet 
(state-'l') occurs at relatively small value of J2/J1 (1.47 
for J3/J1 = 0.3). For 0.5 < J 2 /Ji < (J 2 /Ji) C 2, the 
system is unusual in the sense that the ground state is 
dimer-likc or 5 = 1-like while the singlet-triplet gap is 
filled with many low-lying singlets. For larger values of 
J2, the tctrahcdral chain is equivalent to a standard 5 = 
1 chain not only for the ground state but also for low- 
lying excitations. 

We summarize the results obtained in this section in 
Fig.|, I and Table I. In FigJ| we can see clearly the col- 
lapse of a huge number of singlet states mentioned above. 
Among them, the longest one (which is not shown here) 
comes down from infinitely high energies and hits the RD 
ground state at J2/J1 = {J2/ J\)ci- In Fig.||, some triplet 
branches appear discontinuous at the transition. This 
is because the corresponding state acquires an extensive 
dispersion when the critical value J2 = J-2c is crossed 
from the RD side. The difference between a given type 
of excitation in a 'background' of states '1' and '2' respec- 
tively , i.e. the discontinuity appearing in Fig.^J, vanishes 
as the coupling J 3 decreases. For example, the lfwjel cor- 
responding to the gap A^jpjg t has a discontinuity^!! of the 
order Jf. 



III. MAGNETIZATION PROCESS 

In this section, we investigate the magnetization pro- 
cess of the tetrahedral chain for fields up to the saturation 
field, paying particular attention to the magnetization 
plateaus which appear at m z /m sat = 1/2 and are related 
to two different quantum states. Since the real materials 
have J\ w J2 ~ 40 K, these plateaus may be detected in 
high-field magnetization measurements. 

The appearance of two different types of plateaus at 
TO 2 /m sat = 1/2 is apparent already in the limit of isolated 
tetrahedra (i.e. J3 = 0); for J2 < J\ this plateau occurs 
for magnetic fields J\ < H < J\ + J2 and with the system 
either in state '3' or in state '6' (type I plateau), whereas 
for J2 > J\ the plateau occurs for magnetic fields J2 < 
H < 2J2 with the system in state '9' (type II plateau). 

For a discussion of the magnetization process in the 
presence of interaction ( J3 ^ 0), wc start by considering 
the dimer limit J± J2,3 where the plateau of type I is 
realized. It is well-known that the magnetization process 
in thisJimit can be reduced to that of an effective s = 1/2 
modelEj by regarding a triplet on a dimer bond as an 
upward spin (s z = +5 ) and, correspondingly, a singlet 
as downward spin (s z = —h)- The effective Hamiltonian 
obtained in this way is 
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Weff-1= E [^(^? + l+^i + l) + ^^f + l]+^2 E S I S 1+1 
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(10) 



where J xy and J 2Z are 



j I for symmetric model 
I J3 for generalized model 



'Jzz 



J3 for symmetric model 
J3/2 for generalized model 



(11) 



When J2 is sufficiently larger than J3, crystallization of 
upward spins (or, hardcore bosons) occurs and a 'solid' 
phase . . . 333333 ... or . . . 666666 ... is realized. Note 
that both states spontaneously break link-parity (trans- 
lational symmetry is not broken). While this crystalliza- 
tion is obvious for the symmetric model, the situation is 
slightly subtle for the generalized model since the cou- 
pling between tetrahedra is XY-likc. Fortunately, even 
in this case we can show that either of the two is selected 
by a weak ^-coupling. 

In order to excite this 'crystallized' ground state mag- 
netically, a finite amount of energy (oc Ji + J3) has 
to be spent and this leads to the type-I plateau at 
TO 2 /m sat = 1/2. In particular, it is clear from the ab- 
sence of any kind of kinetic terms that the magnetization 
process for the symmetric model is step-likecj under the 
conditions assumed here (Ji 3> Ji-,Jz). 

As is obvious from Fig. |[ for J\ as J 2 two different 
magnetic particles ('3' or '6' and '9') degenerate and com- 
petition between two different plateau phases (type-I and 
II) occurs. To investigate it, we can use a similar method 
to that used in section I for J 2 ~ J\ and J3 <C J\,Ji- 
Let us fix the magnetization m z /m sat = 1/2. Then the 
half magnetized state is dominated by the three types 
of tetrahedra triplets '3', '6' and '9' and we can con- 
struct an effective S = 1 Hamiltonian using these nearly 
degenerate states. We identify the tetrahedron states 
^= G'3') ± |'6')) and |'9') with S z = ±1 and S z = 0, 
respectively, and obtain, in lowest non-trivial order, the 
following effective Hamiltonian: 



+ (j 2 - ^i)E(^) 2 - 



3 + 1' 



(12) 



It is clear that the above model has two phases: when 
J3 is much larger than J 2 — Ji(> 0), the ground state is 



given by a product 



V J V2 



1 (l'3%-± 



; 6')j) while the trivial 



product £3>j|'9')j becomes the ground state when (S z ) is 
dominant. 

If we notice that the above model is in fact classical, it 
is not difficult to know what kind of transition occurs be- 



tween two phases. Although the Hamiltonian ([12]) looks 
complicated, the fact that [(SJ) 2 - (S*) 2 , (S z f] = en- 
ables us to rewrite it in terms of classical spins Oj which 
take three values — 1, 0, and 1. That is, if we identify 
>(|1) ± I — 1)) and |0) with a = ±1 and a = respec- 
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tively and perform the replacement 
a, = {S*f - (Sff , a] 
the Hamiltonian ( |l2] ) reduces to 



J: 



-f E a 3 a 3 + l + ( J 2 - J l) E ° 2 3 



(13) 



(14) 



This is nothing but tho,Blume-Capel model whose phase 
diagram is well knowncJ; for sufficiently low temperatures 
(zero, in our case), there are a ferromagnetic phase with 
all a taking 1 or —1 and the so-called vacancy phase 
where a = at all sites and these two phases are sep- 
arated by a first-order transition. In the language of 
the original quantum spins, the ferromagnetic state cor- 



responds to ®j|'3')j or 



and the vacancy state 



to <8>j|'9')j. Therefore, the aforementioned two plateau 
phases correspond to the ferromagnetic- and the vacancy 
phases, respectively. 

The above argument is based on the lowest-order ef- 
fective Hamiltonian. If we proceed to the next order, a 
new interaction appears: 
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(7 3 a 3 + l( <7 3 - <J 3 + l) 



(15) 



Fortunately, this type of interactions does not affect the 
two phases mentioned above and we may expect that our 
result will remain qualitatively correct in higher orders 
of perturbation theory. 

We determine the phase boundary by perturbation ex- 
pansion. By comparing the energies of both states ob- 
tained by perturbation expansion (up to second order), 
we obtain the following equation for the phase boundary 



T 7 , 1 j 53Jf 



(16) 
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The leading term J\ — J2 + J3 /4 coincides with that given 
by eq.(fL4|). In Fig.^we show the phase diagram for finite 
magnetization m z /m sat = 1/2. Inside the boundary, we 
have parity breaking states (. . . 33333 ... or . . . 66666 . . . ) 
similar to that found in Ref.|25| and a unique parity- 
symmetric state (. . . 99999 . . . ) outside. 

IV. SUMMARY AND DISCUSSION 

Motivated by recent discovery of tellurates 
CU2TC2O5X5 (X=C1 or Br), we considered low-lying 
excitations and the magnetization process of a system of 
coupled tetrahedra built of four spins S = 1/2. 

To investigate possible low-lying singlets and transi- 
tions between several ground states, we derived a classi- 
cal Ising chain in a magnetic field as an effective Hamil- 
tonian. Using it, we showed the growth of binding energy 
between singlet particles triggers a first-order transition 
between the simple rung-dimer (RD) phase and the spin- 
1 dimer (SID) phase. 

Most of excitations in the RD phase are immobile or 
strongly localized due to the geometry (the so-called 'or- 
thogonal dimer' structures). On the other hand, in ad- 
dition to immobile singlets, excitations with much larger 
mobility (with dispersion ~ J3) arc also allowed in the 
SID phase. In both phases, there exists a 'window' in the 
vicinity of the point J\ = J2, where the lowest excitation 
is not a triplet, but many singlet excitations populate the 
gap to the lowest triplet. A similar situation iii-known to 
occur for a frustrated ladder (J 2 = J3 = JifB- In cer- 
tain limits, these singlets can be viewed as bound states 
of dimer excitations. It would be worth mentioning here 
that these low-lying singlets in the SID phase are not in- 
cluded in the Hilbert space of the effective spin-1 chain. 

We also determined the phase boundary by two dif- 
ferent methods (effective Hamiltonian and a numerical 
method) and compared the results to obtain satisfac- 
tory agreement. In particular, for a generalized coupling 
J3 3> J4 wc concluded that a weakly-coupled perfect 
(Ji = J 2 ) tetrahedra is in the RD phase. As is sug- 
gested by this, the effect of couplings between tetrahedra 
(in particular, the resulting ground states) is sensitive 
to the detail of the couplings (symmetric or generalized, 
Ji and J2) and detailed information on the couplings is 
crucial in comparison with experiments. 

One way to distinguish between two phases (RD and 
SID) experimentally would be to use optical probes. For 
example, the Raman operator which creates singlet ex- 
citations belonging to A-representation of S4 is given 
essentially by Hi with J3 = J4. Since it annihilates 
the RD ground state, the scattering intensity should be 
very weak for (cc)-polarization (c-axis is parallel to the 
chain axis), while it creates singlet combinations of two 



'7' tetrahedra in the SID phase. (The Raman oper- 
ator corresponding to B-representation takes the form 
J2i ^(Si,i-S 3ii )-(S 2 ,i+i -S 4 , i+ i), which excites elemen- 
tary singlets and two-triplet bound states both in the 
RD- and the SID phase.) 

In section III, we investigated the magnetization pro- 
cess. There appears a 1/2-plateau in the magnetization 
curve. In the dimer region J\ J2tH, the plateau 
is accompanied by discrete symmetry (parity) breaking 
and is attributed to the ordered state: . . . 3333333 ... or 
. . . 66666666 . . . , whereas in the 'S = V region a parity- 
symmetric state . . . 9999999 . . . realizes. Note that the 
translational symmetry is not broken at all. The transi- 
tion between these two types of plateaus is described by 
a simple pseudo-spin (S = 1) Hamiltonian equivalent to 
the classical Blume-Capel model. 

Finally, we give some brief comments on the ef- 
fects of possible three-dimensional couplings which recent 
experimental suggest are relatively large. We carried out 
preliminary calculation assuming the simplest interchain 
coupling (Jj_) compatible with the crystal structure and 
found the following: (i) the ground state problem is again 
described by the classical (3D) Ising model at least for 
small couplings and (ii) the phase boundary between RD- 
and SID phase is relatively insensitive to the interchain 
coupling Ji. Therefore we may expect that if we are in 
the SID phase when J± = then so are we even for small 
but finite J±(^- 0). However, new phenomena show up 
when the dynamics in the triplet sector is considered. An 
analyses analogous to the one presented in section II B 
suggests that the triplet ('9') gap gets reduced substan- 
tially while the singlet-singlet gap slightly increases as we 
increase J3 and Jl; spin gaps in the effective spin-1 sector 
finally close and three-dimensional antiferromagnetic or- 
dering may take place in the effective spin-1 system (not 
in the original spin-1/2 system). A similar phenomenon 
has been used as a theoretical tij ick i n the context of the 
so-called composite-spin modclsE3Eil. The region of the 
above ordered phase blows up from the transition line 
between SID and Haldane (see Fig.^) and grows as J± is 
increased. 
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(a) 
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FIG. 1: Spin tetrahedron and two models considered in the 
text: symmetric model (a) and generalized model (b). 




FIG. 2: Energy levels (in unit of Ji) of a single tetrahedron: 
singlet (solid), triplet (broken), and quintet (dot-dashed). En- 
ergy is measured from the ground state. Note that the ground 
state changes at Ji = Ji. 
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FIG. 3: Eigenstates of a single tetrahedron. Singlets: 1 (SD) 
and 2 (SID), triplets: a = (3,4,5), f3 = (6,7,8), and 7 = 
(9, 10, 11), quintet: (12,13,14,15,16). Arrows denote singlets 
and dashed lines and ovals triplets. In the dimer picture, '2', 
'9'-'H', and '12'-'16' can be viewed as 2-triplet bound states. 
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FIG. 4: Ground-state phase diagram for the symmetric model 
(J3 = J4). Note that the phase diagram is symmetric un- 
der J2 <-» J3. The phase boundary between RD- and Hal- 
dane phases as obtained from the effective Hamiltonian Wising 
(eq.^j) is shown by a full line and that from numerical data 
(for 16 sites, or N = 8; Ref. [H^) by open squares. The inset 
shows the same data using the parametrization adopted in 
Ref. H 
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FIG. 6: Several low-lying triplets (in unit Ji) obtained by 
perturbation expansion in J3 (some of them are exact) as a 
function of J2/J1. Parameters are the same as in Fig.^. 1- 
triplet band and 2-triplet continuum are shown by light- and 
dark gray regions. Note that the dispersion of the lowest 
triplet suddenly changes at the transition since the ground 
states on both sides are completely different. 
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FIG. 5: Several low-lying singlets (in unit Ji) obtained by 
perturbation expansion in J3 (some of them are exact) as a 
function of J2 (J3 is fixed). We chose J3 = 0.3Ji (hence 
first-order transition occurs at J2/J1 = jc = 0.96904...). 
Elementary singlet and bound states are shown by solid line, 
while scattering states are shown by broken lines. Note that 
we show only a part of the entire singlet spectrum and actually 
a gap between the singlet ground state and the lowest triplet 
(shown by gray lines) is filled up with singlets composed of 
the elementary singlet discussed in section II. A. 1. 
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TABLE I: Energy, spin, degeneracy, and symmetry classifi- 
cation (in terms of S4) of typical excitations. TV denotes the 
number of tetrahedra. 
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FIG. 7: Phase diagram for m z /m sa t = 1/2 obtained by low- 
order perturbation. Inside the boundary, (3,6)-plateau oc- 
curs. Transition points determined by a similar method to 
that used in section I are plotted by small dots (numerical 
data were taken from Ref.M) . A large dot on the symmetric 
line, which separates (3,6pplateau- and non- pla teau phases 
(shown by a thick line), was taken from Rrfrid. From the 
results for the bond-alternating 5 = 1 chainESEj, we believe 
that the non-plateau phase is realized only on the symmetric 
line. Poor convergence around the symmetry axis J2 = J3 
may be attributed to the proximity to the criticality. 



